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ABSTRACT 

Sticking  probabilities  for  oxygen  on  columbium  and  dilute  columbium- zirconium 
alloys  were  derived  from  ratios  of  oxygen  absorption  rates  to  their  corresponding  oxygen 
fluxes,,  Oxygen  test  pressures  were  varied  from  10  to  10  torr  over  temperatures 
between  1158  and  1373  K.  The  absorption  rates  were  derived  from  a  pressure-drop 
technique  and  fluxes  from  constant  pressure  measurements.  The  kinetics  of  the  reaction 
are  discussed  as  well  as  characteristics  of  the  test  method. 
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ABSORPTION  RATE  STICKING  PROBABILITIES  FOR  OXYGEN  ON  COLUMBIUM 
AND  DILUTE  COLUMBIUM-ZIRCONIUM  ALLOYS 

r 

by  Charles  A.  Barrett 
Lewis  Research  Center 

SUMMARY 

q> 

Reaction  rates  in  oxygen  of  columbium  and  dilute  columbium-zirconium  alloys  (up 

to  0„94  percent  Zr)  were  determined  at  temperatures  between  1158  and  1373  K  at  times 

-7-5 

to  236  hours  and  oxygen  pressures  between  10  and  10  torr.  The  ratio  of  these  rates 
to  their  corresponding  oxygen  fluxes  is  the  sticking  probability.  These  rates  were  de¬ 
termined  by  measurement  of  a  pressure  drop  across  a  molecular  flow  conductance  up¬ 
stream  of  the  pumping  specimen.  The  reaction  rates  were  linear  with  time  and  flux  and 
essentially  independent  of  temperature.  This  indicated  that  surface  adsorption  was  rate 
controlling.!  The  rates  increased  with  increasing  zirconium  content  and  were  consider¬ 
ably  lower  when  the  test-chamber  wall  temperature  and  specimen  temperature  were 
similar  (hot -wall  case),  than  when  the  test -chamber  wall  was  much  cooler  than  the 
specimen  (cold -wall  case). 

The  sticking  probabilities  derived  from  a  regression  analysis  of  the  data  were  0. 047 
and  0.  093  for  columbium  and  columbium  -  1-percent  zirconium  in  the  hot-wall  case  and 
0.  077  and  0. 15  for  columbium  and  columbium  -  1-percent  zirconium  in  the  cold-wall 
case,  i 

A  multiple,  linear  regression  analysis  of  the  data  shows  that  the  linear  absorption 
rate  for  dilute  columbium-zirconium  alloys  can  be  expressed  as  log  R-  =  -1.  3326 
+  1.000  log  flux  +  0.  3028  wt  %  Zr  +  0. 2163  X4  where  Rj^  is  the  linear  absorption  rate, 
flux  is  the  oxygen  flux,  Zr  the  weight  percent  zirconium  present  in  the  columbium,  and 
is  the  factor  which  indicates  whether  the  specimen  is  surrounded  by  a  cold  wall  (1) 
or  a  hot  wall  (0).  The  equation  explains  over  95  percent  of  the  total  variability  with  a 
standard  error  of  estimate  of  0. 1414.  This  equation  shows  (1)  a  linear  dependence  on 
flux  (pressure),  (2)  no  detectable  activation  energy,  (3)  a  monotonic  dependence  on 
zirconium  content,  and  (4)  a  dependence  on  whether  X4  =  1  for  cold-wall  conditions  or 
X4  =  0  for  hot -wall  conditions .  j|  — ^7  (p< 

The  sticking  probabilities  and  estimating  equation  are  valid  until  the  oxygen  solu¬ 
bility  is  exceeded  at  the  metal  surface  and  columbium  monoxide  is  formed.  The  time  to 
reach  this  level  can  be  predicted  for  any  given  flux  and  specimen  thickness  if  the  oxygen 
solubility  limit,  oxygen  diffusion  rate,  and  oxygen  sticking  probabilities  are  known. 


INTRODUCTION 


Columbium  and  its  alloys  can  take  into  solution  l^rge  quantities  of  oxygen  and  other 
interstitial  impurities  at  elevated  temperatures.  If  present  in  sufficiently  large  concen¬ 
trations,  these  dissolved  interstitials  tend  to  degrade  the  desirable  properties  of  the 
metal  such  as  resistance  to  corrosive  attack  by  alkali  metals  and  creep  strength  at 
elevated  temperatures,  which  make  such  alloys  of  interest  for  components  of  advanced 
space-power  systems.  This  type  of  contamination  would  be  negligible  in  space,  but  is  an 
important  consideration  in  ground  testing  of  components  constructed  of  such  alloys. 

These  would  be  proof  tested  in  laboratory  vacuum  environments  where  contaminating 
molecules  are  always  present. 

If  it  is  assumed  that  oxygen  is  the  most  likely  serious  contaminating  gas,  then  the 
maximum  pickup  rate  by  the  alloy  as  a  limiting  case  is  just  the  number  of  oxygen 
molecules  colliding  with  the  surface  (i.e. ,  flux).  For  example,  calculations  based  on  the 
Langmuir  equation  (ref.  1)  indicate  that  an  alloy  tube  with  a  wall  thickness  of  0. 038 
centimeters  would  reach  a  level  of  2000  ppm  oxygen  in  9.62  hours  at  10“  torr  and 
96.2  hours  for  10"6  torr.  This  is  for  the  gas  hitting  the  outer  surface  only.  This  cal¬ 
culation  assumes  that  all  the  gas  in  the  system  is  oxygen  and  that  every  molecule  that 
hits  the  surface  dissociates  and  dissolves  in  the  metal,  oxygen  diffusion  being  fast 
enough  so  that  surface  flux  is  rate  controlling.  However,  if  the  actual  partial  pressure 
of  oxygen  is  used  in  the  flux  calculation  and  if  only  a  fraction  of  the  oxygen  molecules 
colliding  with  the  tube  surface  are  adsorbed,  the  actual  time  to  reach  any  critical  oxygen 
level  would  be  proportionally  longer  or  the  vacuum  requirement  to  reach  the  critical  level 
would  be  that  much  lessened.  The  ratio  of  reacting  molecules  to  colliding  molecules  is 
termed  the  sticking  probability  and  is  presumed  to  be  a  material  constant  for  any  alloy. 

The  purpose  of  this  study  was  to  determine  sticking  probabilities  at  temperatures 
up  to  1373  K  and  oxygen  pressures  of  from  10  ^  to  10  ^  torr  for  pure  columbium  and 
the  nominal  columbium  -  1-percent  zirconium  alloy.  Sticking  probabilities  were  deter¬ 
mined  from  oxygen  pickup  rates  run  at  constant  oxygen  fluxes  at  temperatures  ranging 

from  1158  to  1373  K  for  times  up  to  240  hours. 

In  this  pressure  range  the  anticipated  pickup  rates  were  so  small  that  conventional 
gravimetric  techniques  could  not  be  used  because  they  lack  the  necessary  capacity 
and/or  sensitivity  for  the  total  specimen  weight  change  involved.  Therefore,  the  rate  of 
pickup  was  determined  by  measurement  of  the  pressure  drop  through  a  conductance 
downstream  of  an  oxygen  leak,  with  the  pressure  over  the  specimen,  which  acts  as  a 
pump,  being  constant. 

Other  investigators  who  have  studied  the  low-pressure  oxygen  reactions  of  colum¬ 
bium  or  columbium  -  1 -percent  zirconium  near  1373  K  have  run  at  pressures  greater 
than  ljO-5  torr.  Here,  absorption  takes  place  for  only  a  short  time  before  surface 


scaling  begins  (refs.  2  to  8).  Sticking  probabilities  can  be  derived  from  the  early  part  of 
the  rate  curves  but  are  difficult  to  verify.  Pasternak  has  derived  sticking  probabilities 
for  oxygen  on  pure  columbium  over  the  pressure  and  temperature  ranges  of  interest 
using  a  chemisorption  approach  (refs.  9  and  10).  The  rate  data  from  these  studies  will 
be  compared  with  those  from  this  investigation  in  this  report. 

The  tests  were  performed  at  the  General  Electric  Vacuum  Laboratory,  Schenectady, 
New  York,  and  are  summarized  in  references  11  and  12.  The  data  and  specimens  were 
returned  to  the  NASA  Lewis  Research  Center  for  analyses  and  constitute  the  basis  for 
this  report. 

The  results  of  the  initial  phase  of  this  study  were  published  in  1963  (ref.  13);  since 
then,  the  test  procedure  and  analyses  have  been  revised  to  incorporate  several  improve¬ 
ments  . 


EXPERIMENTAL  METHOD,  APPARATUS,  AND  PROCEDURE 

The  sticking  probabilities  determined  in  this  study  are  simply  the  ratios  of  two  other 
quantities  derived  from  experimental  measurements,  namely,  the  oxygen  absorption  rate 
and  the  corresponding  oxygen  flux  incident  on  the  specimen.  The  techniques  used  for 
determination  of  these  two  quantities  are  described  in  the  following  sections,  along  with 
a  discussion  of  the  experimental  uncertainties  associated  with  each  of  them. 


Oxygen  Absorption -Rate  Measurements 

In  this  study  oxygen  absorption  rates  of  the  test  alloys  were  determined  by  the 
pressure-drop  technique.  This  method  was  originally  developed  by  Wagener  for  the 
molecular  flow  regime  (ref.  14)  in  1950  and  extended  by  della  Porta  (refs.  15  and  16)  and 
Orr  (ref.  17).  The  sensitivity  of  this  technique  can  be  of  the  order  of  10"**  gram,  much 
beyond  the  sensitivity  of  commonly  used  gravimetric  techniques  which,  at  best,  are 
close  to  10”  gram  (ref.  1).  This  method  is  based  on  the  pumping  action  produced  when 
certain  gases  are  brought  into  contact  with  getter  materials.  The  flow  of  gas  being 
gettered  is  measured  by  the  magnitude  of  the  pressure  drop  across  a  known  vacuum  con¬ 
ductance.  The  volumetric  gettering  or  pumping  rate  Q-  can  be  expressed  by  the  follow¬ 
ing  equation: 


Qi  =  apc0  (1) 

where  AP  is  the  measured  pressure  drop  through  an  orifice,  capillary,  or  tube  length 
in  the  molecular  flow  range  of  conductance  CQ.  The  specific  rate  Ri  is 
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Ri=  — 

1  A 

where  K  is  a  conversion  factor  and  A  is  the  specimen  area.  (More  complete  defini¬ 
tions  of  the  symbols  are  presented  in  appendix  A. )  In  this  study  the  rate  was  de¬ 
termined  as  a  function  of  temperature,  pressure  and  time. 

In  order  to  be  able  to  use  the  pressure-drop  method,  it  is  necessary  to  assure  that 
the  measured  pressure  drop  is  attributable  only  to  the  pumping  action  of  the  specimen. 
Experimentation  verified  that  this  was  the  situation  in  our  experiments;  when  no  heated 
specimen  was  present  in  the  system,  there  was  no  detectible  pressure-drop  across  the 
conductance  CQ. 


Apparatus 

The  apparatus  for  studying  continuous  oxygen  absorption  of  columbium  alloys  at  low 
steady-state  pressures  at  temperatures  near  1375  K  was  developed  at  the  General 
Electric  Engineering  Laboratory,  Schenectady,  New  York,  specifically  for  this  program. 
Detailed  summaries  of  the  apparatus  and  experimental  procedure  are  given  in  refer¬ 
ences  10  and  11. 

A  photograph  of  a  typical  test  set-up  is  shown  in  figure  1,  and  a  schematic  of  the 
essential  features  of  the  system  is  shown  in  figure  2.  The  apparatus  consists  of  an 
orifice  throttled,  oil-diffusion-pumped  system,  an  oxygen  permeation  leak  between  the 
pump  and  the  test  section,  a  section  of  known  molecular  flow  conductance,  a  heated  test 
specimen  section,  appropriately  located  vacuum  gages,  and  a  mass  spectrometer.  The 
system  is  constructed  entirely  of  borosilicate  glass  with  a  quartz  test  section.  The  use 
of  an  all  glass  system  was  necessary  in  order  to  avoid  extraneous  pumping  by  metal 
surfaces  other  than  the  specimen.  It  did,  however,  introduce  the  possibility  of  silicon 
transfer  to  the  test  specimen,  which  would  invalidate  the  measurements.  This  possibility 
was  checked  visually  and  by  mass  spectrographic  analysis  of  test  specimens  for  silicon. 
No  evidence  of  silicon  transfer  to  the  test  specimens  was  found.  All  runs  except  one 
were  performed  with  the  test  section  in  a  resistance-wound  tube  furnace.  The  remain¬ 
ing  run  was  induction  heated. 

The  rate  of  oxygen  flow  into  the  test  rig  was  adjusted  by  varying  the  temperature  of 
an  oxygen  permeable  silver  tube.  Generally,  once  the  required  temperature  was 
reached,  the  oxygen  leak  rate  was  quite  stable  and  little  adjustment  was  needed  through¬ 
out  a  given  test.  The  temperature  of  the  test  section  was  monitored  with  a  calibrated 
optical  pyrometer. 
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Figure  1.  -  Low-pressure  oxidationrig  (scale  approximately  1/4). 
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Figure  2.  -  Oxidation  flow  conductance  test  rig. 


Pressure  measurements  were  made  with  miniature  Bayard-Alpert  thorium  dioxide - 
iridium  filament  ionization  gages  with  large  tubulations  and  were  recorded  continuously. 
The  ionization  gages  were  run  with  reduced  emission  currents  in  order  to  minimize  the 
pumping  action  of  the  gages.  It  was  determined  experimentally  that  the  pumping  action 
of  the  vacuum  gages  and  mass  spectrometer  were  negligible  compared  with  the  pumping 
action  of  the  specimen  in  oxygen.  The  pressures  were  measured  at  three  stations,  Pq, 
Pj,  and  Pg  (fig.  2).  For  the  purpose  of  this  experiment,  it  is  vital  that  the  pressure 
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drop  (Pq  -  Pj)  be  accurately  measured.  Therefore,  the  individual  ionization  gages  < 

along  with  their  control  panels  were  calibrated  in  oxygen  over  the  test  pressure  range  of  j 

from  10“ 8  to  10"4  torr  against  two  secondary-gage  standards  and  were  periodically  re¬ 
calibrated  for  drift  (refs.  11  and  12).  The  two  secondary  gages,  along  with  their  indiv¬ 
idual  control  panels,  were  calibrated  for  oxygen  in  an  absolute  system.  The  overall 
uncertainty  in  the  oxygen  pressure  estimate  for  each  gage  was  near  5  percent. 

The  species  of  gas  present  in  the  system  were  checked  with  a  calibrated  magnetic-  j 

deflection  mass  spectrometer  when  steady-state  conditions  were  established  and  at  I 

frequent  intervals  thereafter.  In  the  earlier  runs,  the  mass  spectrometer  pickup  was 
located  near  Pj.  Later,  it  was  located  at  P3  to  give  a  more  accurate  measure  of  the 
gases  over  the  specimen.  In  order  to  measure  oxygen  sticking  probabilities  accurately , 
it  is  necessary  that  interferences  from  other  species  not  be  significant.  At  very  low  < 

pressures,  residual  gases  such  as  carbon  monoxide,  water  vapor,  and  hydrogen  may 
constitute  an  unacceptably  large  fraction  of  the  gases  present.  Therefore,  for  deter- 
ruinations  of  oxygen  sticking  probabilities,  it  was  arbitrarily  decided  to  discard  any  set 
of  data  for  which  the  mass  spectrometer  data  indicated  an  oxygen  concentration  less  than 
85  percent  over  the  specimen.  In  this  concentration  range  (100  to  85  percent)  the  overall 
scatter  in  the  data  was  not  a  function  of  the  oxygen  percentage  (see  appendix  B). 

j 

\ 

j 

Test  Specimens  j 

i 

\ 

A  description  of  the  individual  test  specimens  is  given  in  table  I.  They  were  all  | 

sized  to  have  a  total  surface  area  of  approximately  20  square  centimeters  and  were  cut 
from  sheets  of  several  thicknesses  that  ranged  in  size  from  0. 076  to  0. 127  centimeter. 

The  as-received  sheet  were  j 

(1)  Vapor  degreased  in  trichloroethylene 

(2)  Sonically  cleaned  in  detergent  i 

(3)  Rinsed  in  tap  water  and  then  in  distilled  water 

(4)  Lightly  etched  in  55-percent  water  -  25-percent  nitric  acid  -  20-percent  hydro¬ 

fluoric  acid  (3  to  5  sec)  ^ 

(5)  Rinsed  in  distilled  water  and  air  dried  in  an  oven  at  43  C  (110  F) 

Early  specimens  were  1.27  wide  by  7. 62  centimeters  long;  later  specimens  were  ap¬ 
proximately  2.  54  wide  by  3.  81  or  4.44  centimeters  long.  The  shorter,  wider  specimen 
tended  to  give  a  smaller  pressure  gradient  along  its  length  and  rested  nearer  the  center 
plane  of  the  test  section  than  the  longer  narrower  specimen.  The  derived  rate  data  from 
either  size  did  not  show  any  detectable  difference,  however.  The  surface  finish  of  each 
specimen  is  listed  in  table  I. 

The  tests  were  performed  in  three  different  test  rigs  which  are  described  in  refer- 
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TABLE  I.  -  TEST  SPECIMEN 


[Nominal  alloy  composition,  columbium  -  1  wt  %  zirconium.  ] 


(a)  As -received 


Specimen 

Thickness, 

cm 

Initial  chemistry,  ppm  by  weight 

Surface 

finish, 

jucm 

Sample 

Oxygen 

Nitrogen 

Carbon 

Hydrogen 

Columbium  -  0.75  wt%  zirconium 

930 

180 

20 

4 

152.4 

50,36,37,88 

Columbium  -  0.89  wt  %  zirconium 

360 

23 

48 

6 

106.7 

90 

Columbium  -  0.94  wt  %  zirconium 

30 

64 

27 

1 

38.1 

59,61 

Pure  columbium 

.076 

210 

41 

16 

6.7 

30.5 

52,51 

Pure  columbium 

.114 

92 

55 

32 

3 

53.3 

3 

(b)  Treated 


Sample 

Dimensions, 

cm 

Surface  area, 

cm 

Final  weight, 

g 

Final  oxygen 
content, 
ppm  by  weight 

Number  of 

analyses 

Conductance, 

liter/sec 

3 

1.298  by  7.640 

21.88 

9.56671 

15  900 

11 

0.615 

36 

1.275  by  7.625 

21.77 

10.72576 

3  040 

10 

.615 

37 

1.290  by  7.625 

21.93 

10.  51443 

3  250 

10 

.615 

50 

2.062  by  4.463 

20.10 

10.16088 

2  220 

8 

.615 

51 

2.050  by  4.458 

19.26 

5.58592 

254 

7 

.536 

52 

2.050  by  4.448 

19.20 

5.61180 

5  077 

6 

.615 

59 

2.672  by  3.815 

21.35 

6.53260 

3  600 

10 

.595 

61 

2.672  by  3.815 

21.37 

6.60003 

7  690 

12 

.646 

88 

2.436  by  3.802 

20.11 

10.17256 

5  520 

6 

.619 

90 

2.451  by  3.810 

20.11 

9.03159 

3  330 

7 

.538 

ences  10  and  11.  Their  major  differences  were  in  the  diameter  of  the  test  section  and 
in  the  size  and  position  of  the  conductance  tube.  The  value  for  this  conductance  CQ  at 
its  approximate  temperature  during  each  test  is  also  listed  in  table  I. 


Procedure 

A  test  specimen,  previously  cleaned  and  weighed  to  the  nearest  0.01  milligram, 

was  placed  in  the  apparatus  at  the  end  of  the  test  section.  The  system  was  sealed, 

evacuated,  and  outgassed  at  723  K  overnight.  The  quartz  test  section  was  then  heated  to 

—  ft 

the  test  temperature  and  further  outgassed  until  a  pressure  in  the  10"  torr  region  or 
lower  was  obtained.  The  specimen  was  magnetically  pushed  into  the  center  of  the  test 
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TABLE  n.  -  EXPERIMENTAL  CONDITIONS  AND  TEST  RESULTS  FOR  OXYGEN  ABSORPTION  STUDIES 


Sticking 

probability, 

€ 

i  i  i 

i  i  i 

i  i  i 

1  ■ 

i  i  i 

0.0675 

.0475 

.0651 

.0810 

.0925 

.0525 

0.0525 

.  163 

.0890 

.0783 

.  102 

.108 

0.0701 

.0793 

.0957 

.0846 

0.0275 

0.0415 

.0441 

.0729 

.0737 

.0300 

.0393 

.0511 

.0624 

Oxygen 

flux, 

F, 

(g/cm2)/hr 

i  i  i 
i  i  i 
i  i  i 
i  i  i 
i  i  i 
i  i  i 
i  i  i 
i  i  i 
i  i  i 

in  iji  in  in  to 

i  i  i  i  i  i 

o  o  ©  o  o  o 

1— 1  1— I  1— 1  1— 1  1— I  1— 1 

X  X  X  X  X  X 

E-  CO  03  co  © 

CM  i-l  O  ©  03  © 

CM  CM  CM  H  CD  © 

in  ©  in  ^  if  1? 

i  i  i  i  i  i 

©  ©  ©  ©  o  © 

1— 1  i-H  i— t  i-H  i-t  i— 1 

©  CM  CM  O  CO  Tf 

H  C0  if  rl  H  M 

©  ©  Tf  Tf 

1  1  1  1 
o  o  ©  © 

rH  i—l  i-H  rH 

g  g  g  g 

o  CO  CO  © 

CM  CM  rH  rH 

© 

© 

i—l 

g 

© 

cm’ 

©©©Tf©©TfTf 

1  1  1  1  1  1  1  1 
oo©o©o©o 
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section  where  it  rested  in  a  horizontal  position  with  its  edges  supported  by  the  quartz 
test  section  walls.  The  specimen  was  outgassed  until  the  system  pressure  dropped  close 
to  what  it  was  prior  to  the  introduction  of  the  specimen.  This  sometimes  required  up 
to  three  days  of  pumping.  At  this  time  the  oxygen  leak  was  activated. 

After  the  oxygen  leak  was  activated,  it  typically  took  a  few  minutes  for  P^  to  rise 
to  lxl 0-6  torr  and  for  a  constant  pressure  difference,  for  (PQ  -  Pj),  to  be  obtained.  Zero 
time  was  taken  when  P ^  reached  the  designated  test  pressure. 

At  the  end  of  the  test  period,  the  power  to  the  silver  leak  heater  was  shut  off  as 
well  as  the  power  to  the  specimen  heater.  Shortly  thereafter  P^  dropped  to  below 
lxlO-7  torr  and  no  significant  pressure  difference  was  observed  between  Pq,  Pj, 
and  Pg. 

After  the  apparatus  had  cooled  to  room  temperature,  it  was  backfilled  with  nitrogen, 
it  was  opened,  and  the  specimen  was  removed.  The  specimen  was  reweighed,  sectioned, 
and  analyzed  for  oxygen  by  the  vacuum  fusion  technique.  Selected  sections  of  the 
specimen  were  mounted  for  metallography  and  microhardness  determinations.  In  cases 
where  there  was  apparent  surface  oxidation,  X-ray  diffraction  techniques  were  used  to 
identify  the  surface  oxides . 

In  some  of  the  tests  specimens  (3,  59,  and  51)  the  temperature  and  pressure  condi¬ 
tions  were  maintained  constant  throughout  the  entire  run.  In  the  other  tests,  designated 
as  step  runs,  pressure  and  temperature  conditions  were  successively  varied  stepwise, 
with  the  same  specimen.  On  each  step  change  the  conditions  were  held  constant  long 
enough  to  allow  steady-state  measurements  to  be  taken. 

Test  times,  specimen  dimensions,  surface  area,  average  temperature,  and  average 
pressures  over  the  test  period  of  the  specimen  are  summarized  in  tables  I  and  II. 

Besides  test  times  and  temperature,  the  observed  data  include  test  pressures  at 
Pq,  Pj,  and  Pg.  These  were  recorded  continuously  and  were  arbitrarily  listed  for 
every  hour.  These  values  were  then  converted  to  oxygen  pressures  from  the  calibration 
values.  The  averages  of  these  converted  values  are  listed  in  table  II  for  a  given  test 
temperature  and  time  interval. 

Also  listed  in  table  n,  wherever  possible,  is  the  percent  oxygen  based  on  at  least 
one  mass  spectrometer  scan.  The  location  of  the  mass  spectrometer  and  the  number 
scans  that  were  averaged  is  also  indicated. 


RESULTS 

Determination  of  Oxygen  Absorption  Rates  and  Estimated  Accuracy  of  Determinations 
The  gettering  rate  at  any  time  can  be  calculated  from  the  relation 
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<P0  -  P1>C0K 


(2) 


Here,  Pq  and  are  oxygen  pressures,  CQ  is  the  conductance  between  the  two  pres¬ 
sures,  K  is  a  conversion  factor  based  on  the  ideal  gas  law  to  give  grams  per  hour,  and 
A  is  the  area  of  the  test  specimen.  The  average  rate  for  a  given  time  interval  can  be 
found  by  substituting  average  pressures  in  equation  (2).  The  average  value  for  each  test 
interval  Ft  is  listed  in  table  II.  This  value  was  actually  determined  for  each  hour  of 
the  test  interval  in  order  to  derive  the  standard  deviation  as  well  as  the  average  value. 
The  standard  deviation  divided  by  the  average  value  gives  the  coefficient  of  variation 
mR-,  which  is  a  measure  of  the  variability  of  over  the  test  interval. 

During  any  given  step  or  run,  once  test  conditions  were  set  they  held  quite  steady. 
This  is  indicated  by  the  low  coefficient  of  variation  on  the  rate  constants  which  overall 
averaged  just  under  7  percent,  based  on  the  hourly  values. 

It  is  clear  from  equation  (2)  that  the  accuracy  of  the  oxygen  absorption  rate  deter¬ 
minations  is  strongly  dependent  on  the  accuracy  with  which  the  conductance  CQ  between 
the  gages  at  PQ  and  Pj  is  known  and  on  the  accuracy  of  measurement  of  the  pressure 
drop  between  these  two  stations  (PQ  -  Pj),  as  well  as  on  an  accurate  knowledge  of  the 
specimen  area.  For  the  two  later  test  setups,  used  for  the  greater  part  of  this  study, 
the  conductance  could  be  calculated  in  a  straightforward  way  (ref.  18)  and  is  believed 
to  be  known  to  an  accuracy  of  about  1  percent.  In  the  first  rig  (ref.  11),  the  conductance 
was  experimentally  calibrated  due  to  its  irregular  geometry  and  is  believed  to  be  accu¬ 
rate  to  about  1.  5  percent. 

The  accuracy  of  the  pressure  drop  term  (PQ  -  Pj)  depends  on  how  good  and  stable 
the  oxygen  sensitivity  corrections  are  for  each  gage.  The  uncertainty  associated  with 
each  oxygen  pressure  measurement  was  estimated  as  being  close  to  5  percent. 

The  specimen  area  term  in  the  denominator  of  equation  (2)  is  simply  the  geometric 
area  of  the  specimens.  From  conventional  dimensional  measurements,  this  term  is 
known  with  an  uncertainty  of  about  0.4  percent.  From  these  uncertainties  an  error 
analysis  gave  an  overall  uncertainty  of  close  to  7  percent  for  the  oxygen  absorption  rate. 


Accumulative  Weight  Gains 


From  the  successive  hourly  values  within  a  test  interval,  the  accumulative 
weight  gain  is: 

i=n-l 


W 


p/A 


y~^Ri+Ri+i 

i=0 


At 


(3) 
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where  n  is  the  last  incremental  hour  time  point  and  At  is  the  time  interval.  The  Wp 
values  are  listed  in  table  II.  If  Wp  values  for  time  intervals  of  1  hour  are  plotted 
against  time  on  a  log-log  plot  for  each  test  run  (or  step),  typical  accumulated  weight  - 
gain  curves  like  those  found  in  the  oxidation  literature  result.  Several  such  curves  are 
shown  in  figure  3. 


Time,  hr 

Figure  3.  -  Oxygen  weight  gain  for  various  columbium  and  dilute  columbium-zirconium  alloys  as  function  of  time, 
average  oxygen  flux,  and  specimen  temperature.  Weight  is  nearly  linear  with  time  and  increases  with  increasing 
flux  for  a  given  alloy. 
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Experimental  Verification  of  the  Accuracy  of  Calculated  Oxygen  Absorption  Rates 


As  described  in  the  preceding  section,  the  accumulative  weight  gain  of  a  gettering 
specimen  exposed  to  oxygen  at  high  temperatures  over  a  prolonged  time  period  may  be 
predicted  from  a  knowledge  of  the  oxygen  absorption  rates  over  suitable  time  intervals. 
The  accuracy  of  the  prediction  and  thus  of  the  method  involved  may  be  experimentally 
verified  by  determining  the  weight  change  resulting  from  the  high -temperature  exposure, 
assuming  that  the  total  weight  change  is  due  to  absorption  of  oxygen.  The  validity  of 
this  assumption  can  be  ascertained  by  chemical  analysis  of  the  specimen. 

Such  a  comparison  of  the  predicted  accumulative  weight  gains,  based  on  determina¬ 
tions  of  the  oxygen  absorption  rates  as  described  previously,  with  the  weight  gain  deter¬ 
mined  by  weighing  the  specimen  before  and  after  test  and  finally  with  the  amount  of 
oxygen  increase  during  the  test  as  determined  by  vacuum  fusion  analysis  is  shown  in 
table  IE.  Analysis  of  variance  of  the  data  in  table  m  showed  no  significant  difference  in 
the  values  derived  by  the  three  different  methods.  Over  the  range  of  total  weight  gains 
studied  the  values  agree  to  within  5  percent.  This  agreement  is  considered  excellent 
and  verifies  the  accuracy  of  the  experimental  method  for  the  determination  of  oxygen 
absorption  rates  at  low  oxygen  pressures.  (This  method  can  also  be  used  in  the  con¬ 
trolled  doping  of  oxygen  (or  other  reactive  gases)  for  selected  metal  alloys  at  elevated 
temperatures.  See  ref.  19.) 

There  was  some  uncertainty  as  to  whether  the  long  period  of  outgassing  at  the  test 
temperature  just  prior  to  starting  the  oxygen  leak  changed  the  initial  specimen  weight 
and  oxygen  chemistry  significantly  which,  in  turn,  would  affect  the  AW^.  and  AWC 


TABLE  m.  -  OVERALL  WEIGHT  GAIN  ESTIMATES 
FOR  TEST  SPECIMENS 


Specimen 

Weight  gain,  g 

By  weight  change, 

AW, 

g 

By  oxygen  analysis, 
AWC 

By  pressure  drop, 

AWP 

3 

0.15299 

0.15133 

0.15375 

36 

.01918 

.02264 

.01659 

37 

. 02372 

.02438 

.02069 

50 

.01710 

.01304 

.01454 

51 

.00071 

.00024 

.00035 

52 

.02633 

.02731 

.02350 

59 

.02104 

.02203 

61 

.04381 

. 04924 

.05031 

88 

. 04768 

. 04664 

'  .05385 

90 

.02820 

.02676 

.02435 
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values  listed  in  table  III.  In  order  to  evaluate  this  possibility,  an  0. 127  centimeter  thick 

columbium  columbium  -  0.  75 -percent  zirconium  sample  was  outgassed  for  20  hours 

after  the  usual  baking  out  of  the  system  and  then  removed  from  the  rig  prior  to  testing. 

-4  -8 

The  pressure  dropped  from  10  torr  initially  to  the  mid  10  torr  range  after  20  hours. 
There  was  no  detectable  change  in  weight  or  oxygen  chemistry  of  the  specimen. 


Oxygen  Flux  Over  Specimen 

In  the  preceding  sections,  it  was  shown  that  the  rate  of  oxygen  absorption  at  low 
pressures  can  be  measured  with  high  accuracy  for  gettering  materials  like  columbium 
and  columbium  -  1 -percent  zirconium  by  use  of  the  pressure -drop  technique.  In  order 
to  determine  sticking  probabilities,  it  is  also  necessary  to  know  the  corresponding  oxy¬ 
gen  flux  incident  on  the  specimen.  This  quantity  was  derived  from  the  free  molecular 
flow  expression  (ref.  1) 


F  =  209. 88  P 


(4) 


where  F  is  the  oxygen  flux,  T2^  is  the  absolute  temperature  characteristic  of  the 
oxygen  molecules  in  the  test  region,  M  is  the  molecular  weight  of  oxygen,  and  P2  is 
an  estimate  of  the  oxygen  pressure  over  the  specimen.  Flux  values  for  the  experiments 
previously  described  are  shown  in  table  n. 

Unfortunately,  for  the  experimental  arrangement  used  in  this  type  of  experiment  the 
flux  cannot  be  determined  with  the  same  high  certainty  as  the  rate  of  oxygen  absorption. 
This  stems  primarily  from  the  uncertainty  regarding  the  value  to  assign  to  P2-  An 
estimate  of  the  temperature  of  the  molecules  T2^  is  also  difficult  to  estimate,  but 
fortuitously,  regardless  of  its  value  it  cancels  out  because  it  is  also  used  in  the  deriva¬ 
tion  of  P2-  This  estimate  of  P2  is  derived  from  both  P^  and  P^  and  is  discussed  in 
detail  in  appendix  B.  The  uncertainty  in  this  estimate  of  P2  was  on  the  order  of 
12  percent.  The  uncertainty  in  the  flux  measurement  is  difficult  to  estimate  directly.  If 
equation  (4)  is  assumed  and  an  exponential  drop  in  pressure  along  the  specimen  length 
is  reasonable,  with  the  temperature  of  the  molecules  essentially  cancelled  out,  then  a 
12 -percent  uncertainty  can  be  expected. 


Determination  of  Sticking  Probabilities 

The  ratio  of  the  oxygen  absorption  rate  R.  to  its  corresponding  flux  (average  or 
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instantaneous)  can  be  thought  of  as  the  pumping  efficiency  or,  in  terms  of  equilibrium 
chemisorption  theory,  as  the  steady-state  sticking  probability  (i.e. ,  the  fraction  of 
impinging  molecules  that  actually  stick  to  and  react  with  the  metal  surface).  The  un¬ 
certainty  in  this  estimate,  based  on  the  uncertainties  in  Rj  and  F  was  calculated  to  be 
about  13  percent  derived  from  the  individual  values  of  the  absorption  rates  and  fluxes 
shown  in  table  II . 

Sticking  probabilities  derived  from  multiple  regression  analysis  of  the  data  have 
considerable  utility  because  they  can  be  estimated  for  various  levels  of  the  significant 
variables  along  with  their  upper  and  lower  standard  deviation  limits.  This  method  also 
indicates  which  factors  were  not  significant.  Table  IV  lists  these  derived  sticking 
probabilities  at  selected  levels  of  the  significant  variables  along  with  their  error  limits. 
These  values  were  constant  over  the  flux  range  and  independent  of  temperature  over  the 
range  studied  (1158  to  1373  K).  They  are  strictly  a  function  of  the  zirconium  content  in 
the  columbium  and  whether  the  tests  were  hot  walled  or  cold  walled.  They  range  from 
approximately  0.  05  for  pure  columbium  to  0.  09  for  the  columbium  -  1-percent -zirconium 
alloy  for  the  hot -wall  case.  For  the  cold-wall  case  these  values  are  increased  to  0.  08 
for  pure  columbium  and  0. 15  for  the  columbium  -  1 -percent -zirconium  alloy. 


TABLE  IV.  -  STICKING  PROBABILITIES  DERIVED  FROM  REGRESSING  EQUATIONS  WITH 
UPPER  AND  LOWER  ONE  SIGMA  VALUES 


[Number  of  data  points  included  in  regression  analysis,  46;  =  1.000.] 


Temperature 

effect 

Hot-wall  case,  =  0 

Cold-wall  case,  =  1 

Standard 

deviation 

Columbium 

Columbium  -  1-percent 

zirconium 

Columbium 

Columbium  -  1 -per cent 

zirconium 

Temperature,  K 

1158 

1373 

1158 

1373 

1158 

1373 

1158 

1373 

Sticking  probabilities 

Included 

0.0617 

.0444 

.0319 

0.0656 

.0472 

.0340 

0.125 

.0900 

.0648 

0.133 

.0957 

.0689 

0.104 

.0747 

.0537 

0.210 

.151 

.109 

0.224 

.161 

.116 

+1  sigma 

-1  sigma 

Excluded 

0.0644 

.0465 

.0336 

0.129 

.0934 

.0674 

0.106 

.0765 

.0552 

0.212 

.154 

.111 

+1  sigma 

-1  sigma 
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Regression  Analysis  of  Rate  Data 


Separate  log-log  plots  were  made  for  pure  columbium  and  the  nominal  columbium  - 
1-percent -zirconium  alloy  of  the  linear  rate  constants  against  their  respective  fluxes 
for  each  test  condition.  These  are  shown  in  figures  4  and  5  for  the  furnace -heated 
specimens.  Specimen  88  was  induction  heated;  data  for  this  run  are  plotted  separately 
on  figure  6.  The  temperature  of  each  test  is  indicated  also.  Although  there  is  a  fairly 
broad  scatter  band,  the  data  points  seem  to  fall  on  nearly  straight  lines  with  a  slope  of 
near  unity.  The  maximum  possible  value  for  sticking  probabilities  is  unity  (i.  e. ,  Ri 
=  flux).  A  line  representing  this  line  is  also  plotted  in  each  of  the  figures. 

A  simple  linear  regression  was  run  for  the  data  on  each  plot  following  the  model 

log  Ri  =  a  +  bj  log  flux  (5) 


Figure  4.  -  Linear  oxygen  absorption  of  pure  columbium  as  function  of  oxygen  flux 
and  specimen  temperature  for  hot-wall  case.  Diagonal  lines  indicate  range  of 
sticking  probabilities  from  regression  equation.  The  regression  line  is  derived 
from  nine  data  points. 
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Figure  5.  -  Linear  oxygen  absorption  of  dilute  columbium-zirconium  for  hot-wall 
case.  Diagonal  lines  indicate  range  of  sticking  probabilities  from  regression 
equation.  The  regression  line  is  derived  from  30  data  points. 


TABLE  V.  -  LINEAR  REGRESSION  VALUES:  Y  =  a  +  b1X1  +  b9X9  +  b,,X«  +  b4X 


Standard  deviation  of  regression  response  surface. 
'Dropped,  insignificant  to  a  probability  level  of  0.  95. 


Because  the  value  of  in  each  case  was  close  to  unity  (table  V),  the  equation  was  re¬ 
calculated  with  this  slope  forced  to  unity.  The  error  was  increased  only  slightly  as 
measured  by  the  standard  error  of  estimate.  The  fraction  of  explained  variation  is  also 
high  in  all  the  three  cases,  which  verifies  the  validity  of  equation  (5)  as  a  good  first 
approximation. 

If  all  the  data  for  the  three  plots  (46  data  points)  are  combined  and  run  in  the  same 
model  (eq.  (4)),  the  fit  is  poorer,  and  the  error  correspondingly  higher  (table  V). 


Multiple  Linear  Regression  Analysis  to  Isolate  Effects  of  Flux,  Temperature, 
Zirconium  Content  and/or  Other  Variables 

Because  there  was  a  larger  degree  of  unexplained  variation  when  all  the  data  points 
were  combined  in  equation  (5),  it  seemed  reasonable  to  assume  that  there  were  other 
factors  that  significantly  affected  the  total  variation.  Temperature  and  zirconium  con¬ 
centration  seemed  the  most  obvious,  although  it  appeared  the  cold-wall  induction  heated 
rates  were  higher  than  the  corresponding  tests  run  in  the  hot -wall  furnace.  These  other 
factors  along  with  flux  could  be  evaluated  by  Multiple  Linear  Regression  Analysis 
(ref.  20)  with  a  suitable  computer  program.  (The  program  used  for  this  study  was  the 
Multiple  Linear  Regression  FORTRAN  IV  Program,  RAPIER  written  by  S.  Sidik, 

B.  Henry,  and  R.  Lambird  of  Lewis). 

With  this  technique  an  equation  for  a  response  surface  can  be  estimated  which  un¬ 
mixes  the  various  factors.  This  occurs  because  the  coefficients  are  derived  from 
partial  differential  equations  where  the  partial  derivative  of  the  total  variation  is  taken 
with  respect  to  one  factor  while  all  others  are  held  constant,  hence,  the  term  partial 
regression  coefficients.  In  so  doing,  the  degree  of  fit,  scatter,  and  the  partial  regres¬ 
sion  coefficients  for  various  groupings  of  the  rate  data  can  be  analyzed. 

The  following  model  equation  was  used: 

y  =  a  +  bjXj  +  b2X2  +  bgXg  +  b4X4  (6) 


where 


y  =  log  Rj 
Xj  =  log  flux 
X2  =  1/T2 
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X3  =  wt  %  Zr 


and 


X4  =  0  or  1 

This  equation  represents  the  effects  of  flux,  temperature,  percent  zirconium,  and 
a  fourth  factor  X4,  which  is  a  dummy  variable  used  to  represent  whether  the  tests  were 
induction  heated  (X4  =  +1)  or  otherwise  (X4  =  0).  The  last  factor,  if  it  is  significant, 
merely  shows  how  much  higher  the  rate  is  in  the  cold-wall  case.  When  the  data  are 
broken  into  smaller  groups  just  representing  pure  columbium  and  one  type  of  heating, 
then  the  equation  reduces  to  just  the  first  two  variables. 

In  the  final  analysis  all  46  rate  constants  could  be  combined  as  a  single  equation. 

The  flux  coefficient  b^  was  numerically  close  to  one  and  was  forced  to  unity  with  only  a 
slight  increase  in  error. 

Table  V  lists  the  constants  for  both  b4  floating  and  bj  forced  to  unity.  The  regression 
program  is  designed  so  that  each  of  the  coefficients  is  given  a  statistical  test  of  signifi¬ 
cance  to  a  given  probability  level.  For  example,  we  can  be  at  least  95  percent  sure  that 
any  factor  we  choose  to  keep  is  a  real,  not  just  a  random  effect.  In  this  study  a  95- 
percent  level  was  used,  and  any  coefficient  testing  less  than  this  was  automatically 
dropped  and  the  equation  was  recalculated.  In  the  46 -point  calculation  temperature  was 
dropped  on  this  basis.  This  insignificance  can  also  be  inferred  qualitatively  if  the  ratio 
of  a  coefficient  to  its  deviation  is  small.  The  final  estimative  equation  was 

log  R.  =  -1.  3326  +  1.000  log  flux  +  0.  3028  percent  Zr  +  0.2163  X4  (7) 

This  equation  explains  over  95  percent  of  the  total  variation  with  a  standard  error  of 
estimate  of  0. 1414. 

The  sticking  probabilities  listed  in  table  IV  were  derived  from  this  final  equation. 

Any  given  flux  can  be  substituted  into  the  equation  to  estimate  its  corresponding  reaction 
rate  R^.  The  standard  error  of  estimate  is  added  or  subtracted  before  the  antilog  is 
taken  to  get  the  upper  and  lower  limits  of  R- .  These  values  divided  by  the  corresponding 
flux  give  the  sticking  probabilities.  These  equations  estimate  most  precisely  over  the 
levels  investigated.  Extrapolations  are  usually  questionable  unless  theoretical  consider¬ 
ations  warrant  them.  For  example,  there  should  be  no  lower  theoretical  limit  on  flux, 
but  extrapolations  to  5  percent  zirconium  would  be  questionable.  Also,  tests  at  higher  or 
lower  temperatures  might  well  give  significantly  different  rates. 


/ 
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Comparison  with  Other  Investigations 

Other  investigators  who  have  run  absorption -rate  studies  on  columbium  or  other 

works  from  which  such  rate  constants  can  be  derived  from  the  early  stages  of  low- 

pressure  oxidation  studies  are  Hurlen,  et  al.  (ref.  2),  Inouye  (ref.  3),  Spruiell  and 

Harms  (ref.  4),  Pasternak  and  Evans  (ref.  10),  Kofstad  and  Espevik  (ref.  5),  and 

Delgrosso,  et  al.  (refs.  6  and  7).  The  columbium  -  1  percent -zirconium  alloy  was 

studied  by  Hogan,  et  al.  (ref.  8).  These  data  cover  a  wide  range  of  oxygen  pressures 
-7-2 

from  10  to  10  torr.  For  our  purposes,  these  pressure  values  were  converted  to 
oxygen  fluxes  using  equation  (5)  after  using  the  transpiration  correction  (eq.  (Bl)  on 
their  pressure  estimates.  The  temperature  at  station  one  Tj  was  assumed  in  all 
cases  to  be  298  K.  Figure  7  shows  the  linear  rate  absorption  constants  plotted  against 
flux  for  all  the  reference  values,  as  well  as  the  data  from  this  investigation. 

Because  so  many  points  were  involved,  the  30  points  from  figure  5  in  this  investi¬ 
gation  were  represented  by  their  unforced  regression  line.  The  data  of  reference  7, 
made  up  of  67  data  points,  were  also  plotted  in  this  manner.  Here,  two  lines  were  used 
to  represent  the  two  extreme  temperatures  (1143  and  1478  K)  because  temperature  did 
significantly  change  the  rates. 

The  temperatures  indicated  on  the  plot  varied  from  1123  to  1478  K.  Delgrosso  and 
Hogan  ran  all  their  tests  with  an  induction  heater  so  these  can  be  considered  cold-wall 
tests.  Pasternak  used  a  columbium  ribbon  heated  by  self-resistance  which  also  can  be 
considered  a  cold-wall  test.  The  remainder  were  hot -wall  tests,  except  for  our  speci¬ 
men  88. 

Considering  the  wide  ranges  of  temperature  and  pressure  covered,  most  of  the  data 
seem  to  agree  fairly  well,  particularly  for  pure  columbium.  There  does  appear  to  be  a 

change  of  slope  for  the  pure  columbium  data  from  near  unity  at  the  low  fluxes  to  near 

-3 

0. 6  at  fluxes  near  10  grams  per  square  centimeter  per  hour,  and  higher.  The  effect 
of  temperature  is  large  in  the  case  of  the  reference  7  data,  perhaps  because  of  the 
larger  temperature  span.  The  most  puzzling  aspect  of  this  plot  are  the  high  values 
Hogan  found  in  his  columbium  -  1 -percent  zirconium  runs.  These  would  give  e’s 
approaching  unity.  His  rates  are  probably  accurate  because  they  were  obtained  with  a 
microbalance.  His  oxygen  pressure  estimates  could  be  low  due  to  the  pumping  action  of 
the  metal  chamber  walls.  In  all  other  cited  data,  the  test  rigs  were  constructed  of 
either  glass,  quartz,  or  ceramic.  A  narrow  tubulation  to  an  ionization  gage  with  a 
significant  pumping  rate  could  also  be  a  factor. 

In  general,  figure  7  shows  that  the  sticking  probabilities  for  the  cold-wall  runs 
were  higher  than  the  hot -wall  runs.  Also  zirconium  additions  gave  higher  sticking 
probabilities.  In  summary,  the  literature  data  show  the  same  general  trends  as  our 
data  except  for  a  more  pronounced  temperature  effect. 
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Figure  7.  -  Linear  oxygen  absorption  rates  for  coiumbium  and  dilute  columbium-zirconium 
alloys  for  various  investigations  as  function  of  estimated  oxygen  flux  and  specimen  tem¬ 
perature.  {Includes  individual  data  points  and  regression  line  fit  of  many  data  points. 
Numbers  with  each  data  point  are  specimen  temperatures  in  K. 


DISCUSSION 


Mechanism  of  Oxidation 


Characteristics  of  the  low-pressure  oxidation  of  coiumbium  and  the  nominal 
coiumbium  -  1 -percent -zirconium  alloy  at  elevated  temperatures  have  been  described 
by  Kofstad  and  Espevik  (ref.  5),  Barrett  and  Rosenblum  (ref.  13),  Pasternak  (ref.  9), 
Inouye  (ref.  3),  and  others  in  their  low-pressure  oxidation  studies  on  coiumbium.  The 
kinetics  that  they  found  can  be  divided  into  three  main  stages:  (1)  an  initial  linear 
rate  stage,  (2)  a  subsequent  stage  with  a  decreasing  oxidation  rate,  and  (3)  a  final 
stage  characterized  by  an  increasing  rate  of  oxidation. 
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The  initial  linear  stage  is  due  to  the  solution  of  oxygen  after  chemisorption  (i.  e. , 
adsorption)  on  the  surface.  Here,  the  rate  controlling  step  is  the  oxygen  flux  on  the 
metal  surface  which  is  primarily  a  function  of  the  oxygen  pressure,  provided  the  metal 
temperature  is  high  enough  so  that  diffusion  into  the  metal  is  essentially  instantaneous. 

If  the  oxygen  diffusion  rate  is  high  enough,  the  metal  surface  is  effectively  clean  (i.  e. , 
coverage  is  low). 

Stage  2  commences  when  the  diffusion  rate  can  no  longer  keep  up  with  the  adsorption 
rate.  This  occurs  almost  instantaneously  at  high  pressures  unless  the  specimen  is 
quite  thick.  The  solubility  limit  is  then  reached  and  surface  oxidation  (stage  3)  com¬ 
mences. 

In  our  investigation  the  test  pressures  were  low  enough  so  that,  for  the  tempera¬ 
tures,  thicknesses,  and  times  investigated,  the  intitial  linear  rate  stage  was  always 
maintained.  The  rate  curves  for  all  the  steps  in  table  II  were  similar  to  those  in  fig¬ 
ure  3  which  is  typical  of  stage  I  oxidation.  The  linear  rate  constants  derived  for  each 
step  should,  if  surface  adsorption  is  rate  controlling,  be  directly  proportional  to  oxygen 
pressure  or  more  precisely  to  oxygen  flux.  Figures  4,  5,  and  6  show  that,  indeed,  the 
linear  rate  constants  are  directly  proportional  to  oxygen  flux.  These  plots  and  regres¬ 
sion  analyses,  however,  also  show  that  there  is  no  strong  temperature  dependence  for 
our  data,  at  least  between  1158  and  1373  K. 

Pasternak  gives  the  most  plausible  explanation  of  the  independence  of  the  sticking 
probability  on  temperature  (ref.  9).  The  rate  of  absorption  is  strictly  a  function  of  the 
arrival  rate  of  the  oxygen  (i.  e. ,  flux)  and  the  number  of  dual  sites  available  on  the  sur¬ 
face.  The  molecules  are  chemisorbed  and  diffusion  of  atomic  oxygen  follows.  No 
significant  physisorption  takes  place  at  the  test  temperatures  involved.  Activated 
desorption  is  apparently  not  involved. 

It  can  be  inferred  from  calculations  that  the  chemisorbed  oxygen  diffuses  into  the 
metal  much  faster  at  low  pressures  than  it  arrives  both  for  pure  columbium  and 
columbium -zirconium  alloys,  based  on  the  activation  energy  for  diffusion  at  these  rela¬ 
tively  high  test  temperatures.  Thus,  for  these  tests  the  metal  surface  is  always  effec¬ 
tively  clean  (i.  e.  ,  coverage  is  low).  Under  these  conditions,  no  significant  activation 
energy  is  involved  and  the  sticking  probability  is  constant. 

Two  questions  which  have  yet  to  be  resolved  are  (1)  why  small  amounts  of  zirconium 
significantly  raise  the  sticking  probability  when  alloyed  to  columbium  and  (2)  why  cold- 
wall  runs  can  show  up  to  twice  the  sticking  probability  of  the  hot -wall  runs  ?  At  this 
time,  no  explanation  can  be  given  for  the  zirconium  effect;  a  possible  explanation  of  the 
second  question  is  suggested  as  follows: 

Four  different  investigators  ran  cold -Wall  tests  (fig.  7).  They  were  Pasternak  and 
Delgrosso  for  pure  columbium  and  Hogan  and  this  investigation  for  columbium  - 
1 -percent  zirconium  (sample  88  in  this  study).  In  all  three  investigations,  which  used 
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induction  heating,  the  sticking  probabilities  were  much  higher  than  for  the  hot -wall  case. 
In  the  remaining  cold -wall  series  of  tests,  Pasternak  used  self -resistance  heating  of  the 
columbium  and  here  the  sticking  probabilities  were  also  somewhat  higher  than  those 
derived  for  the  hot -wall  case. 

The  reasons  for  these  differences  in  the  hot -wall  and  cold -wall  sticking  probabilities 
are  not  known.  This  whole  experimental  technique  is  an  averaging  process  for  the 
oxygen  molecules  in  the  system.  It  would  seem  that  the  individual  behavior  of  the 
molecules  would  be  critical  as  to  how  the  electronic  levels  of  the  oxygen  atoms  are  in¬ 
volved  in  the  adsorption  process. 

An  observation  made  by  Fleming  (private  communication  to  the  author  by  R.  B. 
Fleming,  General  Electric  Co. ,  Research  and  Development  Center,  Schenectady,  N.Y. , 
October,  1965)  was  that  sticking  probabilities  can  be  expected  to  drop  as  the  kinetic 
energy  (i.  e. ,  temperature)  of  the  colliding  molecules  increases.  He  suggests  that  a 
factor  of  2  or  3  difference  in  kinetic  energy  between  room -temperature  molecules  and 
1273  K  molecules  is  not  unreasonable.  The  exact  relation  would  depend  on  the  exact 
temperature  of  the  furnace  walls  and  specimen,  the  location  and  number  of  molecular 
collisons  before  absorption,  and  the  accommodation  coefficients  for  oxygen  on  hot 
quartz  and  hot  columbium. 


Limitations  of  Use  of  Sticking  Probabilities  -  End  of  Stage  I  Oxidation 


As  long  as  surface  adsorption  is  the  rate  determining  step,  sticking  probabilities 
can  be  used  to  estimate  oxygen  contamination  of  getter  metals  like  columbium.  This 
implies  that  the  rate  of  diffusion  of  oxygen  into  the  metal  is  faster  than  the  effective 
arrival  rate  (e  x  flux)  of  oxygen  to  the  metal  surface.  This  can  easily  be  verified  from 
jump -frequency  calculations  based  on  the  activation  energies  of  the  diffusion  of  oxygen 
into  the  alloy,  (ref.  9). 

This  surface  flux  would  be  expected  to  be  rate  controlling  until  the  oxygen  concen¬ 
tration  at  the  metal  surface  exceeds  the  solubility  limit  for  oxygen  in  columbium  at  the 
temperature  of  interest.  At  this  point  surface  scaling  begins  and  stage  I  oxidation 
(absorption)  ends. 

When  a  constant  surface  flux  process  is  rate  controlling,  the  solution  of  Fick's  dif¬ 
fusion  differential  equations  for  the  proper  boundary  conditions  gives  the  following 
equation  (ref.  21): 


C 


F  Dt  .  3y2  -  l2 
D  ,2  g2 


oo 

L  V  -I"  e-^n2n2t/l2  cog  iwy 

?Lj  n2  * 

n=l 


(8) 
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For  values  of  Dt /f2  >  1 


the  series  term  is  negligible  and  the  equation  reduces  to 


(9) 


where  is  the  initial  concentration,  C  is  the  concentration  at  anytime  t,  Fe  is  the 
effective  (i.  e. ,  flux  X  e)  surface  flux,  D  is  the  diffusion  coefficient,  t  is  the  time,  l  is 
the  distance  from  the  surface  to  the  center  of  the  specimen,  and  y  is  any  point  meas- 
red  from  the  center. 

This  is  the  nonsteady-state  diffusion  case  and  predicts  relatively  shallow  oxygen 

o 

concentration  gradients.  The  gradient  depends  only  on  Tit/ 1  .  For  example,  for  pure 
columbium  at  1253  Kwith  l  =  0.038  centimeter  (0.015  in.)  and  t  =  1  hour,  this  ratio 
equals  1.039  and  gives  a  Cy_z/Cy_0  of  1.53.  For  t  =  10  this  Dt/Z2  ratio  is  10.39, 
which  gives  a  concentration  ratio  of  approximately  1. 04.  In  our  tests  for  the  total 

O 

specimen  exposure  time,  Dt /l  was  always  much  greater  than  one,  indicating  a  shallow 
gradient. 

Microhardness  traverses  across  the  specimens  showed  near  constant  hardness 
through  the  specimens,  which  indicates  a  lack  of  gradient.  Removal  of  surface  layers 
by  grinding  and  oxygen  analysis  of  the  remainder  also  indicated  no  apparent  gradient. 

Under  this  condition  the  surface  concentration  of  oxygen  will  rise  gradually  from  the 
initial  bulk  oxygen  concentration.  The  gradient  will  keep  the  same  shape  as  the  surface 
concentration  rises  until  this  concentration  reaches  C  the  solubility  limit  for  the 
formation  of  columbium  monoxide  at  a  given  temperature.  When  Cy_^  =  Cg,  surface 
oxide  will  start  to  form.  Until  this  time,  the  process  will  be  characterized  by  the  rate 
relations  described  previously:  shallow  oxygen  concentration  gradients  and  no  detectable 
surface  oxides.  Apparently  the  surface  coverage  is  low  even  close  to  C_.  This  is  true 
even  in  the  nominal  columbium  -  1-percent -zirconium  alloys  where  the  initial  saturation 
of  the  zirconium  as  zirconium  monoxide  or  zirconium  dioxide  causes  no  apparent  change 
in  the  oxidation  kinetics.  Columbium  oxide  surface  film  formation  is  again  critical. 

Once  the  surface  oxide  is  sufficiently  thick,  it  is  assumed  the  constant  surface  flux 
will  no  longer  be  rate  controlling.  Here,  the  process  will  become  controlled  by  diffu¬ 
sion  through  the  growing  oxide  and  will  be  proportional  to  the  square  root  of  time  (i.  e.  , 
parabolic). 

If  this  hypothesis  is  correct,  it  should  be  possible,  for  a  given  temperature  and 
thickness  of  material,  to  predict  when  the  change  of  kinetics  from  a  linear  to  a  parabolic 
rate  will  occur.  Equation  (9)  can  be  rewritten  with  time  as  the  dependent  variable  as 
follows: 


(Cs-Ci>  Z2 
Fe  3D 


(10) 
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where  t  is  time  for  saturation,  F  is  the  flux,  and  C  is  the  solubility  limit  for 
S  o 

columbium  monoxide  at  y  =  Z,  Ci  is  the  initial  oxygen  concentration,  and  e  is  the 
sticking  probability  for  pure  columbium.  Beyond  this  tg  limit,  different  sticking 
probabilities  involving  the  surface  oxide  would  apply,  although  they  would  probably 
not  be  rate  controlling. 

The  validity  of  this  t  estimate  was  verified  experimentally  for  a  pure  columbium 

s 

sample  (see  appendix  C). 

Thus,  when  a  getter  alloy  like  columbium  is  exposed  at  elevated  temperatures  to 
oxygen  in  the  molecular  flow  regime,  the  oxygen  sticking  probability  is  the  critical 
material  constant  in  predicting  oxygen  contamination. 


SUMMARY  OF  RESULTS  AND  CONCLUSIONS 

The  investigation  of  absorption  rate  sticking  probabilities  for  oxygen  on  columbium 
and  dilute  columbium -zirconium  alloys  yielded  the  following  conclusions: 

1.  The  derived  low-pressure  steady-state  sticking  probabilities  (ratio  of  absorption 
rate  to  flux)  between  1158  and  1373  K  for  oxygen  on  pure  columbium  and  the  columbium  - 
1  percent  zirconium  alloy  are  0. 047  and  0.  093,  respectively.  These  results  are  for 
conditions  when  the  temperature  of  the  colliding  oxygen  molecules  is  nearly  the  same  as 
for  the  specimen  (i.e. ,  hot -wall  case). 

2.  Under  comparable  test  conditions,  but  with  the  temperature  of  the  colliding  oxygen 
molecules  relatively  low  compared  with  the  1158  to  1373  K  specimen,  the  sticking 
probability  is  0.  077  for  pure  columbium  and  0. 15  for  columbium  -  1  percent  zirconium 
(i.e.,  cold -wall  case). 

3.  The  reaction  rate  data  indicate  that  surface  oxygen  adsorption  is  the  rate  control¬ 
ling  process  giving  linear  absorption  rates  with  time  and  oxygen  flux,  and  no  significant 
activation  energy.  The  linear  rate  constant  increases  with  increasing  zirconium  content 
and  for  a  given  oxygen  flux  tends  to  increase  with  decreasing  kinetic  energy  of  the 
oxygen  molecules. 

4.  A  multiple  linear  regression  analysis  of  the  oxygen  absorption  rate  data  showed 
that  the  following  estimating  equation  can  be  used: 

log  Rj  =  -1.  3326  +  1. 000  log  flux  +  0.  3028  wt  %  Zr  +  0.  2163  X4 

where  Rj  is  the  linear  absorption  rate,  flux  is  the  oxygen  flux,  Zr  is  the  zirconium 
content  in  the  columbium  base  alloy,  and  X4  is  the  factor  which  indicates  whether  the 
specimen  is  surrounded  by  a  cold  wall  (+1)  or  a  hot  wall  (0).  The  equation  explains  over 
95  percent  of  the  total  variability  with  a  standard  error  of  estimate  of  0. 1414. 


27 


5. 1  The  rate  and  sticking  probability  estimates  are  only  valid  when  the  oxygen  sur¬ 
face  flux  is  rate  controlling.  This  is  no  longer  rate  controlling  when  the  solubility  limit 
of  oxygen  in  the  alloy  is  exceeded  on  the  specimen  surface  forming  columbium  monoxide 
The  time  to  reach  this  value  t  can  be  predicted  from 

S 

.  <CS-Ci>  !2 


where  Cg  is  the  solubility  limit  of  oxygen  in  the  alloy  when  columbium  monoxide  is 
formed,  Ci  is  the  initial  oxygen  concentration,  F  is  the  oxygen  flux,  e  is  the  sticking 
probability  of  oxygen  on  the- alloy,  l  is  one-half  the  specimen  thickness,  and  D  is  the 
oxygen  diffusion  coefficients.  f) 

Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Cleveland,  Ohio,  June  13,  1968, 

129-03-03-01-22. 
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APPENDIX  A 
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Li 


l 


M 


mR, 


Qi 

Ri 


Ai 

T2(G) 

T2(s) 


SYMBOLS 

specimen  area,  cm^ 

K 

O 

concentration  of  oxygen  in  alloy  at  any  time  t,  g/cm 

O 

initial  concentration  of  oxygen,  g/cm 
molecular  flow  conductance,  liters/sec 

O 

oxygen  saturation  limit,  g/cm 

oxygen  diffusion  coefficient  in  alloy,  cm  /hr 

inside  diameter  of  test  section,  cm 
n 

oxygen  flux,  (g/cm  )/hr 

effective  oxygen  flux  (e  x  oxygen  flux) 

standard  deviation  of  R^  over  test  interval 

conversion  factor  from  (torr)(liters)/sec  to  g/hr,  =  (M/Tq)  x  57.  729 

distance  from  to  leading  edge  of  specimen,  cm 

lengths  of  i  test  section,  cm 

1/2  thickness  of  specimen,  cm 

molecular  weight  of  reacting  gas,  g 

coefficient  of  variation,  crR^/R^ 

oxygen  pressure  at  ionization  gage  station  0,  torr 

oxygen  pressure  ionization  gage  at  ionization  gage  station  1,  torr 

estimated  average  oxygen  pressure  over  specimen 

oxygen  pressure  at  ionization  gage  station  3,  torr 

volumetric  pumping  rate,  (torr)(liter)/sec 

average  oxygen  pickup  rate  over  test  interval,  (g/cm  ) /hr 

temperature  at  station  0,  K 

temperature  at  station  1,  K 

estimated  temperature  of  the  colliding  oxygen  molecules,  K 
temperature  at  test  specimen,  K 


t 

WP/A 

y 


time,  hr 

2 

accumulative  weight  gain  over  the  test  interval,  g/cm 
total  specimen  weight  gain,  g 
any  distance  along  thickness,  t,  cm 
oxygen  sticking  probability  (rate/flux) 
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APPENDIX  B 


PROBLEMS  ASSOCIATED  WITH  ESTIMATING  OXYGEN  FLUX  OVER  SPECIMEN 


The  flux  over  the  specimen  was  estimated  by  equation  (3)  as  F  =  209.  8  P2  yM/T 2(g)- 
This  depends  on  the  accuracy  of  estimating  P2  and  T2-  Because  it  was  not  feasible  to 
measure  this  pressure  directly  with  a  gage,  it  was  determined  from  (1)  the  downstream 
gage  P3  and  (2)  the  upstream  gage  P^  The  downstream  Pg  oxygen  pressure  cor¬ 
rected  for  temperature  (i.  e. ,  the  transpiration  correction)  is  (ref.  1) 


(B-l) 


where  Pg  is  the  downstream  oxygen  pressure,  T2^  is  the  absolute  temperature 
around  the  specimen,  Tg  is  the  absolute  temperature  at  Pg  and  Pg_g  is  the  estimated 
oxygen  pressure  at  the  downstream  end  of  the  sample.  The  uncertainty  of  this  estimate 
was  0.  05.  The  second  estimate  from  Pj  is  more  complicated  because  it  not  only  in¬ 
volves  a  transpiration  correction  from  the  cool  zone  to  the  furnace  zone,  but  also  a 
pressure  drop  due  to  the  pumping  action  of  the  specimen  and  the  conductance  drop  be¬ 
tween  Pj  and  the  beginning  of  the  specimen.  This  upstream  pressure  is  (ref.  12) 


pl-2  '  P1  - 


„  /Vi^tLYt 


1/2 


2(G) 


QiT2(G) 
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3.81  DjTq 


3.81  T 
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i=2 


where  the  first  term  is  the  pressure  drop  in  the  borosilicate  glass  tube  just  after  Pj 
times  the  transpiration  correction  on  entering  the  hot  section.  The  second  term  repre¬ 
sents  the  additional  pressure  drop  through  the  various  diameter  quartz  sections  to  the 
leading  edge  of  the  specimen.  The  and  are  the  length  and  diameter,  respec¬ 
tively,  of  the  borosilicate  glass  tube.  The  1^  and  Di  are  the  length  and  diameter  of 
the  i^1  quartz  segment,  the  flow  rate  Qj  equals  (Pq  -  Pj)Cq,  and  Tq,  T-^,  and 
are  the  temperatures  in  various  test  sections.  The  total  length  L  from  P.^  to  the 
leading  edge  of  the  specimen  is  also  tabulated  for  each  run.  All  these  values  are  listed  in 
table  VI.  The  overall  uncertainty  of  this  estimate,  assuming  the  validity  of  the  equation, 
is  about  11  percent,  based  on  the  uncertainties  due  to  Pj,  Qj,  and  D.  The  derived 
values  of  Pg  2  and  Pj  2  are  tabulated  in  table  II  for  the  dimensions  and  temperatures 
involved  for  the  various  specimen  runs.  If  the  two  pressures  (Pg_2  and  are 

equal,  then  there  is  no  measurable  pressure  drop  across  the  specimen.  In  this  study 
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TABLE  VI.  -  PARAMETERS  TO  ESTIMATE  Pl2  FOR  OXYGEN  (M  =  32) 


Specimen  Oxygen  pressure  Volumetric 
at  station  1,  pumping  rate 


Absolute  temperature,  K 


—  Station  0,  Station  1,  Oxygen 

^i’  Tq  Tj  molecule, 

(torr)  (liter)  T  /G\ 


Dimensions  of  test  setup,  cm 

First  segment  Second  segment  Third  segment  Length  of 


Pj  2  was  always  significantly  greater  than  Pg  gj  which  implies  a  pressure  drop  along 
the  specimen.  According  to  Fleming  (ref.  11)  this  drop  is  exponential  in  nature.  There¬ 
fore,  to  estimate  the  average  pressure  over  the  specimen,  the  logs  of  Pg_2  and  Pj_2 
are  averaged  and  its  antilog  estimates  P2  are  used  in  equation  (3).  (In  the  earlier  work 
(ref.  13),  there  was  no  Pg  gage,  so  a  pressure  drop  was  not  accounted  for  in  the  flux 
calculation.  Instead,  oxygen  flux  over  the  specimen  was  simply 


Flux  =  209.88 


\l/2 

rH 

T2(G)/ 


These  gave  biased  high  flux  values  and  corresponding  biased  low  sticking  probabilities. 
(Introduction  of  the  Pg  downstream  gage  introduced  the  whole  problem  of  pressure  drop 
to  and  across  the  specimen  as  well  as  indicating  extraneous  gases  over  the  specimen. ) 

Obviously,  the  closer  the  Pj_2  and  Pg_g  are  to  their  estimates,  the  more  con¬ 
fidence  in  the  P2  estimate.  The  ratios  varied  from  1. 09  to  8.  90  for  columbium- 
zirconium  alloys  with  an  average  of  3.  36.  For  pure  columbium  the  ratio  of  Pj_2 
P3  g  varies  from  1. 47  to  2. 05  with  an  average  of  1.60.  This  drop  is  a  function  of  the 
diameter  of  the  test  section,  the  length  of  the  specimen,  and  its  sticking  probability. 

Vacuum  fusion  analysis  for  oxygen  of  0.635  centimeter  long  segments  along  the 
specimen  length  also  indicate  an  oxygen  pressure  gradient  along  the  specimen.  This 
gradient  could  be  inferred  from  these  values,  except  for  oxygen  diffusion  along  the 
specimen  as  well,  which  would  tend  to  flatten  the  gradient. 

These  derived  average  fluxes  are  listed  in  table  II  for  each  interval.  The  estimated 
uncertainty  in  deriving  P2  and,  hence,  the  flux  is  about  13  percent  based  on  the  un¬ 
certainties  in  determining  Pg_2  and  P^ _2-  The  validity  of  the  flux  equation  depends  on 
the  validity  of  the  transpiration  correction  which  is  used  in  both  equations  (Bl)  and  (B2). 

The  transpiration  correction  applies  strictly  to  the  molecular  flow  regime  when  the 
Knudsen  number  Kn  is  quite  large  (>10),  and  a  state  of  equilibrium  exists  as  defined  by 
Knudsen  as  zero  net  mass  transfer  across  a  plane  rather  than  an  equality  of  pressure 
(ref.  22).  In  our  study  the  Knudsen  number  was  quite  high  ranging  from  50  to  500  but, 
because  of  the  pumping  action  of  the  specimen,  Knudsen's  equality  of  flux  is  only  an 
approximation. 

Another  problem  in  using  equation  (3)  is  in  using  the  proper  temperature  T2^. 
This  is  really  a  measure  of  the  temperature  (i.  e. ,  kinetic  energy)  of  the  oxygen  mole¬ 
cules  colliding  with  the  specimen  surface.  The  oxygen  molecules  enter  the  system  at  a 
temperature  T^  and  presumably  collide  at  least  once  with  the  chamber  wall  before 
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striking  the  specimen.  If  the  accommodation  coefficient"^  is  unity,  one  collision  is 
sufficient  in  the  hot -wall  case  to  raise  the  temperature  to  The  cold-wall  case 

is  more  complicated  because  the  molecule  is  bouncing  back  and  forth  between  a  hot 
specimen  and  the  colder  wall,  but  it  always  strikes  the  specimen  with  the  energy  of  the 
wall,  approximately  in  our  case. 

It  is  more  likely  that  the  accommodation  coefficients  are  not  unity,  however,  so  it 
appears  the  oxygen  molecules  striking  the  specimen  surface  are  always  at  some  inter¬ 
mediate  temperature,  probably  much  lower  in  the  cold-wall  case  than  in  the  hot -wall 
case. 

Another  complicating  factor  was  in  estimating  Pg  in  equation  (Bl)  because  of  the 
presence  of  extraneous  gases  detected  at  this  gage  by  the  mass  spectrometer.  These 
gases  were  primarily  carbon  monoxide  and  carbon  dioxide,  but  water  vapor,  helium  and 
hydrogen  were  also  present.  These  gases  were  found  in  increasing  amounts  with  de¬ 
creasing  pressures  (refs.  11  and  12).  This  was  not  so  serious  a  problem  at  Pj,  which 
was  usually  at  least  a  decade  higher  than  Pg.  This  showed  up  as  an  upward  bias  from  a 
straight  line  in  a  P^  against  Pg  log -log  plot.  No  bias  was  apparent  at  a  P^  >  3x10  ^ 
torr.  Values  less  than  this  that  were  unbiased  indicated  at  least  85  percent  oxygen  at 

_7 

the  Pg  gage.  On  this  basis  all  data  with  values  of  Pj  <  3x10  torr  were  discarded 
unless  a  mass  spectrometer  at  Pg  indicated  at  least  85  percent  02-  The  data  points 
in  this  grouping  with  percent  oxygen  estimates  at  the  Pg  gage  were  corrected  by 
multiplying  the  derived  flux  by  the  percentage  of  oxygen  present.  A  regression  analysis 
performed  on  this  corrected  data  indicated  no  significant  difference  in  the  fit  or  in  the 
regression  constants  for  this  data  from  the  same  data  uncorrected  for  oxygen.  Thus, 
for  the  46  data  points  used,  no  percent  oxygen  correction  was  considered  necessary. 
Attempts  to  correct  the  biased  data  for  extraneous  gases  were  not  successful. 

The  source  of  carbon-monoxide  and  carbon  dioxide  is  in  the  reaction  of  oxygen  with 
carbon  in  the  refractory -alloy  specimen  and  gage  filaments.  Helium  diffuses  through  the 
hot  quartz  while  comes  off  the  glass  and  quartz  walls.  Hydrogen  is  outgassed  by 
the  hot  specimen  and  is  a  reaction  product  of  water  with  the  specimen. 


=  (Tg^Q)  -  1’i)/('r2(G)  "  w^ere  a  *s  accommodation  coefficient  ranging 
from  0  to  one,  Tg^  and  T^  are  the  specimen,  and  P^  is  the  temperature  and  Tg^ 
is  the  temperature  of  the  oxygen  molecules  after  colliding  with  the  hot  wall  or  specimen. 
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APPENDIX  C 


EXPERIMENTAL  VERIFICATION  OF  PREDICTED  END  OF  STAGE  I  OXIDATION 


The  time  when  oxygen  absorption  is  no  longer  rate  controlling  can  be  predicted  from 


t 


s 


(Cs-Ci>  _Z2 
F£  3D 


(10) 


where  t  is  the  time  for  surface  saturation  of  oxygen  in  columbium,  l  is  half  the  thick¬ 
ness,  Cg  is  the  solubility  limit  of  oxygen  in  the  metal  gms  at  some  temperature  T,  Cj 
is  the  initial  concentration  of  oxygen  in  the  metal,  D  is  the  diffusion  coefficient  of 
oxygen  in  the  metal  at  T,  and  F£  is  the  effective  oxygen  flux,  gms. 

Because  both  oxygen  diffusion  rates  and  solubility  limits  are  available  for  pure 
columbium  as  a  function  of  temperature,  it  is  possible  to  choose  a  combination  of  test 
temperature  and  oxygen  pressure  where  a  change  from  a  surface -flux  controlled  reaction 
to  a  surface -oxidation  controlled  reaction  should  occur  say  within  a  2  36 -hour  test  run. 

Such  a  test  was  performed  on  a  pure  columbium  sample  0. 114-centimeter  thick 
(Z  x  0. 057  cm)  at  927°  C  (sample  3  in  tables  I  and  n).  The  Ci  was  0. 0092  percent  with  a 
C  calculated  to  be  0.  315  weight  percent  at  927°  C  (ref.  23).  The  value  of  D  used  was 

S  n 

2.  80x10”'  square  centimeter  per  second  (ref.  24)  and  converted  to  square  centimeters 
per  hour.  Because  pure  columbium  has  an  value  of  e  of  0. 047,  it  was  estimated  that  an 
approximate  test  pressure  of  3xl0“5  torr  oxygen  would  cause  Cg  to  be  exceeded  within 
a  236 -hour  test  run,  so  a  break  could  be  detected.  Figure  8  shows  the  accumulated 
weight  gain  plot  against  time  on  a  log -log  plot  of  this  special  rim.  A  change  in  mechanism 
would  show  up  as  change  in  slope  from  one  to  one -half.  It  was  calculated  that  this  break 
should  occur  at  close  to  t„  =  26  hours.  The  actual  P9  test  pressure  was  estimated  to 
be  3. 19x10  torr.  In  this  run  there  was  no  gage  present  at  Pg  so  the  estimated  flux 
was  derived  strictly  from  Pj_2  aad  the  Pi_2/^3-2  rati°s  f°r  the  other  pure  columbium 
tests. 

There  is  an  apparent  break  in  the  curve  somewhere  between  20  and  100  hours  but  it 
is  not  as  sharp  as  expected.  This  is  partly  because,  as  the  rate  started  to  slow  down, 
the  pressure  at  P^  started  to  build  up,  thus,  it  was  difficult  to  hold  Pj  constant.  Dur¬ 
ing  the  test  Pj  rose  from  4.  71  to  6.  00x10”^  torr.  This  compared  with  the  straight  lin¬ 
ear  runs  where  the  Pj  could  be  held  quite  constant  with  a  coefficient  of  variation  less 
than  10  percent. 

The  initial  break  from  a  linear  rate  is  close  to  20  hours.  From  20  to  100  hours 
there  is  an  intermediate  rate  characterized  by  a  patchy  buildup  of  surface  oxide.  The 
terminal  behavior  appears  parabolic. 
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Figure  8.  -  Oxygen  weight  gain  for  pure  columbium  (0. 114  cm  thick)  as  function  of  time  at  1300  K  and  in  10"5  torr 
oxygen  pressure  range.  Curve  shows  weight  gain  to  be  linear  with  time  until  20  hours.  After  100  hours  behavior 
appears  parabolic.  Initial  break  point  can  be  estimated  from  equation  (10). 


When  the  uncertainties  in  diffusion  coefficients,  solubilities,  and  flux  determinations 
are  considered,  it  was  thought  that  the  agreement  with  the  calculated  value  of  26  hours  is 
good.  This  break  to  a  parabolic  rate  was  further  confirmed  by  a  tight  sooty  black  oxide 
completely  covering  the  specimen  showing  that  surface  oxidation  did,  indeed,  take  place. 
The  scale  was  shown  by  X-ray  diffraction  to  be  predominantly  columbium  dioxide  with 
a  trace  of  columbium  pentoxide.  The  oxygen  concentration  gradient  was  also  apparently 
quite  shallow  through  the  specimen  as  indicated  by  the  Knoop  microhardness  traverses. 
This  implies  the  bulk  concentration  was  close  to  Cg  throughout  as  was  expected  when 
process  e  is  initially  rate  controlling. 
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